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ABSTRACT 



Context. Molecular hydrogen (H2) is the most abundant molecule in the circumstellar (CS) environments of young stars, and is a key 
element in giant planet formation. The measurement of the H2 content provides the most direct probe of the total amount of CS gas, 
especially in the inner warm planet-forming regions of the disks. 

Aims. Most Herbig Be stars (HBes) are distant from the Sun and their nature and evolution are still debated. We therefore conducted 
mid-infrared observations of Ht as a tracer of warm gas around HBes known to have gas-rich CS environments. 
Methods. We report a search for the H2 S(l) emission line at 17.0348 ^m in the CS environments of 5 HBes with the high resolution 
spectroscopic mode of VISIR (ESO VLT Imager and Spectrometer for the mid-InfraRed). 

Results. No source shows evidence for H2 emission at 17.0348 fim. Stringent 3cr upper limits on the integrated line fluxes are derived. 
Depending on the adopted temperature, limits on column densities and masses of warm gas are also estimated. These non-detections 
constrain the amount of warm (> 150 A") gas in the immediate CS environments of our target stars to be less than ~ 1 - lOMj^p. 

Key words, stars: circumstellar matter - stars: formation - stars: pre-main sequence - ISM: molecules 



1. Introduction 

Molecular hydrogen (H2) is the most abundant molecule in the 
environment of young stars. It remains optically thin up to high 
column densities (~10^^ cm^), and does not freeze onto dust 
grain surfaces. It furthermore self-shields efficiently against pho- 
todissociation by far-ultraviolet (FUV) photons. Therefore, H2 is 
the only molecule that can directly constrain the mass reservoir 
of molecular gas in the circumstellar (CS) environment of pre- 
main sequence stars. On the other hand, H2 is one of the most 
challenging molecules to detect. Electronic transitions occur in 
the UV to which the Earth's atmosphere is opaque, and rotational 
and rovibrational transitions at infrared (IR) wavelengths are 
faint because of their quadrupolar origin. FUV spectroscopic ob- 
servations of H2 absorption lines have provided evidence for the 
presence of cold and warm excited H2 around numerous Herbig 
Ae/Be stars (HAeBes) but have not allowed us to determine the 
spatia l distribution of the observed gas (e.g. Martin-Zaidi et al. 
l2008h . Detection s of H2 hne em ission from disks around young 
stars with ISO (Thi et al."2001') were contradicted by g round- 
based observations (Richter et a l.]|2002l;[Sako et alj20"05l) . which 
showed that the observed e mission c ould b e due to the surround- 
ing cloud material. Carm ona et alj ( l2008h modeled the mid-IR 
H2 lines originating in a gas-rich disk, seen face-on, surround- 
ing a Herbig Ae (HAe) star at 140 pc from the Sun. By as- 
suming that the gas and dust were well-mixed in the disk, a 
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gas-to-dust ratio of about 100, and that Tgas - Tdusi, those au- 
thors demonstrated that mid-IR H2 lines could not be detected 
with the existing instruments. Indeed, they were unable to de- 
tect any H2 mid-IR emission line in their sample of 6 HAe stars. 
At the present time, observations of H2 mid-IR emission lines 
have been reported in 8 of 76 T Tauri stars observed with Spitzer 
(iLahuis et al J I2007V Although mid-IR windows are strongly af- 
fected by sky and instrument background emission, the advent 
of high spectral and spatial resolution spectrographs allows us 
to study the H2 emission from the ground, as demonstrated in 
the cases of t wo HAeBes HD 97048 (Martin-Zaidi et al. 2007|) 
and AB Aur ( Bitner et al.ll2007l) . for which particular conditions 
are required in their CS disks. The analyses of these data usually 
assumes that the H2 excitation is in local thermodynamic equi- 
librium (LTE) and can thus be characterized by a single excita- 
tion temperature, which should be close to the gas temperature 
because of the low critical densities. The AB Aur observations 
infer an H2 gas temperature that is significantly higher than the 
dust temperature. 

In contrast to HAes, the more distant Herbig Be stars (HBes) 
earlier than B9 type have been poorly studied. As a conse- 
quence, the nature and evolution of the CS material suiTound- 
ing HBes is still a subject of controversy. The Spectral Energy 
Distributions (SEDs) are very different from one HBe to another 
(e.g. Hillenbrand et al. 1992). Some have SEDs that are compa- 
rable to those observed for HAes, and have been classified as 
group I or group II stars haboring disks according to the classi- 
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Table 1. Astrophysical parameters of the sample stars (columns 2 to 6). Viad is the radial velocity of the star in the heliocentric rest 
frame. Columns 7 to 12 summarize the observations. The airmass and seeing intervals are given from the beginning to the end of 
the observations. 



Star 


Sp. 
Type 


(K) 


(mag) 


(kms-') 


d 
(pc) 


(s) 


Airmass 


Optical 
Seeing 

(") 


Standard 
Star 


Airmass 


Optical 
Seeing 
(") 


HD 98922 


B9 


10470<" 


0.34<" 


-15*^' 


>540'" 


2700 


1.14-1.15 


0.76-1.41 


HD 25025 


1.024-1.025 


0.69-0.71 


HD 250550 


B7 


12800<^' 


0.57<3' 


+31W 


606±367<5' 


3600 


1.32-1.50 


0.83-1.50 


HD 93813 


1.018-1.023 


0.91-1.05 


HD 259431 


B5 


15900<3' 


0.88<3' 


+43<''> 


290±84<5> 


2250 


1.23-1.32 


0.75-1.28 


HD 39425 


1.12-1.14 


1.51-1.62 


HD 76534 


B2 


20000<''' 


0.80<" 


+ 17(4) 


>160<" 


3600 


1.05-1.14 


0.85-1.25 


HD 93813 


1.014-1.017 


0.92-0.96 


HD 45677 


B2/B1 


21400*" 


0.87<" 


+21.6<«' 


500'"' 


3450 


1.02-1.08 


0.69-1.40 


HD 25025 


1.024-1.025 


0.69-0.71 


References: van den Ancker et al. ( 1998): Acke et al. (2005): Bouret et al. (2003); Finkenzeller & Jankovics ( 1984j): '^iBrittain et al.l 
(I2007il; *''>iMartin et alj(i2004,): <''>,Valenti et al.,(,200Q): *'*>,Evans, ll967i): C'.Zorec et al.,ll99&) as Quoted in Cidale et aLaOOlb. 



fication of lMeeus et al. I (l200l . However, numerous HBes SEDs 
with little IR excess have been modeled successfully by free-free 
emission originating in a gaseous CS envelope, and compared 
with rapidly rotating classical Be stars. In addition, the mid-IR 
emission observed in HBes is generally not confined to optically 
thick disks but originates in more complex environments such as 
remn a nt envelopes or halos dLeinert et al.ll200ll: iPolomski et alj 
I2OO2I: IVinkovic et"!!] l2006h . This implies that there are struc- 
tural differences between HAes and t he more luminous HBes 
dNatta et al.ll2000l: iLeinertet al]|200ll) . These results are fully 
consistent with the more rapid evolution of the more massive 
HAeBes, which should still be partially embedded in their natal 
cloud. On the other hand, near-IR interferometric d ata of HBes 
have b een interpreted using flared disks models (e.g. lKraus et alj 
l2007h . 

We selected five HBe stars whose CS en vironments are rich 
in gas as observed by the FUSE satellite (iMartin-Zai'di et all 
I2OO8). The selected stars have IR excesses that are assumed to 
repre s ent the presence of CS dusty disks (e.g. lHillenbrand et"aLl 
Il992t IPolomski etani2002h . but the disk scenario has not yet 
been clearly established. The aim of our study is to detect the 
H2 pure rotational emission line at 17.0348 //m in the vicinity 
of these stars with the high resolu tion mode of the VLT/VISIR 
spectrograph dLagage et al.l l2004l) and constrain the warm CS 
gaseous component. 

2. Selection of the target stars 

The selection of the five target stars was motivated by the pre- 
vious observations of H2 electronic lines in the FUV spectral 
range, which showed that the CS environ ments of HBes were 
rich in cold and warm H2 dMartin-Zaidi et al. 2008). The detec- 
tion of H2 absorption lines provided evidence for the presence of 
cold (-100 K) and warm/hot CS gas (up to 1500 K) around these 
stars. We first briefly describe the CS environments of each of the 
five sample stars. The main astrophysical parameters of the stars 
are summarized in Table [1] 

HD 98922: This s tar is classified as a se lf-shadowed disk 
on the basis its SED dvan Boekel et al.ll2003l) . Its optical spec- 
trum displays broad, but relatively weak [O I] emission pro- 
files, which cannot originate in the dusty disk surface but pos- 
sibly in a rotating gjiseous disk inside the dust-sublimation ra- 
dius dAcke et al.li2005ft . Our preliminary analysis of its FUSE 
spectrum indicates that the CS environment of this star is rich 
in H2. The excitation of H2 appears t o be similar to the HBes 
observed by Mart in-Zaidi et al.l d2008l) . The modeling of the ex- 
citation conditions with the Meudon FDR Code (Le Petit et al. 
I2OO6I) allows us to conclude that the observed gas is probably 
located in the remnant cloud where the star formed. We empha- 



size that no H2 mid-IR line was detected by t he Spitzer Space 
Telescope towards this star (iLahuis et al.ll200'7l) . 

HD 250550 and HD 259431: The SEDs of these two 
stars show signs of strong near-IR excesse s that were mod- 
eled w ith flat optically thick accretion disks bv lHillenbrand et akl 
dl992l) . Spectro-polarimetric observations pr ovided indication s 
of flattened structures surrounding the stars dVink et al.ll2002l) . 
However, spectroscopic observations provided no evidence for 
Keplerian rotation with sin gle-peaked emiss ion lines detected 
for both H I Bry and CO (.Brittain et al.ll2007h . The excitation of 
H2 in the CS environments of these stars, revealed by the FUSE 
satellite, favored an interpretation in terms of CS envelopes 
which were the remnants of the parent molecular clouds in which 
the s tars were formed dBouret et al.l 120031 : iMartin-Zaidi et al.l 
I2008h . 

HD 76534: iHiUenbrand eTaTI d 19921) analyzed the SED of 
HD 76534 and stressed that the low near-IR excess was proba- 
bly due to free-free emission in an ionized envelope rather than 
CS dust. By analyzing FUV absorption lines of H2 and deriving 
the excitation conditions of the gas around the star. ' Martin et al.l 
(2004) concluded that the CS environment of HD 76534, which 
is a B2 star, is probably too hostile for the presence of a CS disk. 

HD 45677: The measured infrared excess was attributed t o 
a dust shell surrounding HD 45677 bv lSwings& Allenldl971h . 
The presence of a disk was inferred by polarization and U V spec- 
troscopic measurements of accreting gas (Schulte-Ladbeck et al.l 
ll992l:lGradv etanil993l) : these studies also agreed with the pres- 
ence of an actively accreting CS disk around the star. The pres- 
ence of such a disk is c onfirmed by the analysis of emission lines 
from ionized metals by Muratorio et al.l d2006l) . Our preliminary 
analysis of the FUSE spectrum of this star shows the presence 
of two or more gaseous components along the line of sight. This 
may be consistent with the presence of a disk associated with a 
CS envelope. 



3. Observations and data reduction 

Searches for mid-IR H2 lines can be undertaken with the 
ESO/VISIR instrument due to its high spectral and spatial reso- 
lution. Its high spectral resolution (10000 < R < 30000) allows 
us to disentangle the H2 lines from the absorption lines due to 
the Earth's atmosphere. VISIR off'ers access to the most intense 
pure rotational fines of H2: the S(l) line at 17.0348 yum, S(2) at 
12.2786 yum, S(3) at 9.6649 yum, and S(4) at 8.0250 jum. The S(0) 
transition close to 28 fim is not observable from the ground due 
to the Earth's atmospheric absorption, and the S(3) line, which 
is located amidst a forest of telluric ozone features, is only ob- 
servable for particularly favorable Doppler shifts. 
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Table 2. 3cr upper confidence limits on the integrated fluxes and intensities of the S(l) line, and upper Umits on the total column 
densities of H2 and masses of H2 as a function of the adopted temperature. Agbs is the expected position of the line in the observed 
spectra. 



Star 




Integrated flux 


Intensity 
(ergs cm"- s"' sr"') 


A'(H2) upper limits 




H2 mass upper 


limits"" 


HD 


(jum) 


(ergs s"' cm"^) 




(cm-2) 






"(Mj„p~10-3 


Me) 










150K 


300K 


lOOOK 


150 K 


300 K 


1000 K 


98922 


17.0328 


<1.2xI0-" 


<1.2xl0-^ 


1.1x10-3 


6.9x10^' 


1.6x10-' 


3.3 


2.0x10^' 


5.9x10-^ 


250550 


17.0373 


<1.6xI0-'-* 


<3. 8x10-3 


1.5x10^3 


9.3x10^' 


2.2x10-' 




3.4+5-^x10-' 


1.0+i«xl0-' 


259431 


17.0377 


<1.8xI0-''* 


<4.2xl0-3 


1.6x10^3 


1.0x10^2 


2.4x10-' 


8.6+J]xl0-2 


2.5+|>10-2 


45677 


17.0365 


<1.6xI0-''* 


<3.7xl0-3 


1.4x10-3 


9.1x10^' 


2.1x10-' 


3.7 


2.3x10-' 


6.7x10-2 



(a) Masses of H2 are calculated assuming the distances quoted in Table[lJ(see text). For HD 98922, the lower limit on the distance is used and 
the mass scales with the distance as cf. 



The five sample stars were observed with the high spectral 
resolution long-slit mode of VISIR on January 15 and 16, 2006. 
The central wavelength of the observations was set to be 17.035 
fim. We used the 0.75" slit, providing a spectral resolution of 
about 14000. The details of the observation conditions are sum- 
marized in Table [1] According to its IRAS fluxes, HD 76534 
should be observable in the spectral range of VISIR. However, 
during the observation, the mid-IR flux of this star appeared to 
be very low and the resulting spectrum could not be used to com- 
plete an accurate analysis. 

For all other observations, the standard "chopping and nod- 
ding" technique was used to suppress the large sky and telescope 
background dominating at mid-IR wavelengths (for details of the 
observation technique see Marti n-Zaidi et al.l 120071) . The chop- 
ping was performed with a 8" amplitude, which enabled us to 
detect the gas of any elongated structure around our targets. To 
correct the spectrum for the Earth atmospheric absorption and 
to obtain the absolute flux calibration, we observed standard 
starsQ immediately before and after observing our sources. Since 
the HD 45677 and HD 98922 spectra have far higher signal-to- 
noise (S/N) ratios than that of the standard stars, we divided the 
spectra of HD 98922, HD 250550, and HD 259431 by that of 
HD 45677, and divided the spectrum of HD 45677 by that of 
HD 98922, to correct for the telluric absorption. We checked 
that the Doppler shifts were sufficiently different from one star 
to the other, to ensure that the analysis would not be signifi- 
cantly affected. We used the o bserved and modeled spectra of 
the standard stars (ICohen et al.llT999) to derive the absolute flux 
calibration. The wavelength calibration was completed by fitting 
the observed sky background features with a model of Paranal's 
atmospheric emission. The spectra were not connected for dust 
extinction, which is negligible in the mid-IR for Ay< 40 mag 
(iFluks et al.lll99l . 

4. Data analysis 

For the four high S/N spectra of stars in our sample, we dis- 
play in Fig. [1] the spectra of the standard stars overplotted on 
the target spectra (top panel), the entire atmosphere-corrected 
spectra before the flux calibration (central panel), and the re- 
gion of the flux-calibrated atmosphere-corrected spectra where 
the S(l) H2 line should appear if present (bottom panel). None 
of the observed sources showed evidence for H2 emission close 
to 17.035yum. 

In the flux-calibrated spectra, we used two different meth- 
ods to determine the standard deviation (cr) of the continuum 
flux. First, we calculated the standard deviation for wavelength 
ranges relatively unaffected by telluric absorption, and close to 
the wavelength of interest. Second, we estimated cr in regions of 

' http://www.eso.org/sci/facilities/paranal/instruments/visir/tools/ 



the spectra with comparable atmospheric transmission to the line 
of interest. The first method provided more conservative values 
of cr for all but one star, namely HD 98922; for this star, the line 
of interest should be detected at wavelengths corresponding to 
the strong telluric absorption line, and the atmospheric transmis- 
sion played an important role. For our analysis, we kept the more 
conservative values, and used the cr values obtained with the first 
method for HD 250550, HD 259431, and HD 45677, and the 
value given by the second method for HD 98922. The 3cr upper 
limits on the integrated line fluxes were calculated by integrating 
over a Gaussian of full width at half maximum (FWHM) equal 
to a spectral resolution element (Av ~ 21 kms-'), and an ampli- 
tude of about 3xcrx (flux), centered on the expected wavelength 
for the S(l) line (Fig.[T]). The upper limits on integrated fluxes 
and intensities of the S(l) line derived for each star are tabulated 
in Table |2] These limits were typ ically a factor of 2 lowe r than 
the detection towards HD 97048 dMartin-Zaidi et al.ll2007 ). The 
intensities were calculated by dividing the integrated fluxes by 
the solid angle of the slit of about 0.75" width with a spatial 
resolution element of about 0.427". 

From the upper limits on integrated intensities and by assum- 
ing that the emitting H2 is optically thin at LTE, we estimated 
the total column densities of H2 (for details on the method, see 
van Dishoeck 1992) as a function of adopted temperatures (150, 
300, 1000 K). Even though the FUSE data indicated the presence 
of a gas component with T ~100 K for our target stars, we mea- 
sured 150 K to be the lowest temperature of the gas, which can 
emit in the H2 S(l) line. For temperatures below 150 K, the S(0) 
line at 28 yum (not observable from the ground) provides more 
reliable constraints on the optically thin gas properties than the 
S(l)line. 

Under the same assumptions and assuming a homogeneous 
medium, the mass of warm H2 is given by: 



Fui ^ 
{hc/4nA)A,ix,iT) °' 



where F„/ is the line flux, d is the distance in pc to the star, A 
is the wavelength of the transition u - I, A„i is the spontaneous 
transition probability, x„(r) is the fractional population of the 
level M at the temperature T in LTE, and / is the conversion fac- 
tor required for deriving the total gas mass from the H2-ortho 
or H2-para mass, assuming that ortho /para ratio is con trolled by 
the gas temperature (from Eq.(l) in lTakahashill200l1) . Our re- 
sults are presented in Table |2l The upper limits on the masses 
are calculated by assuming the distances quoted in Table [1] For 
HD 98922, we used the lower limit on the distance and its mass 
scales with the distance as cP. 
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Fig. 1. Spectra obtained for the H2 S(l) line at 17.0348 fim. Top panel: continuum spectra of the standard star and of the target 
before telluric correction. Central panel: full corrected spectra: dotted lines show spectral regions strongly affected by telluric 
features. Bottom panel: zoom of the region where the S(l) line should be observed (dashed vertical lines). A Gaussian of width 
FWHM = 21 km s ' and integrated line flux equal to the 3cr line-flux upper limits is overplotted. The spectra were corrected neither 
for the radial velocity of the targets nor the Earth's rotation velocity. 



5. Discussion 

We observed five HBe stars with the high resolution spectro- 
scopic mode of VISIR to search for H2 pure rotational S(l) emis- 
sion at 17.0348 jim. As probed by previous FUSE observations, 
the CS environments of these stars provide evidence of large 
reservoirs (A^(H2) ~ 10^' cm^) of cold H2 (T ~ 100 K), and also 
evidence for the presence of warm/hot excited H2 {T > 500 K). 
The H2, observed by FUSE, probably corresponds to the rem- 
nant of the molecular cloud in which the stars were formed 
and the observed H2 FUV lines do not originate in a CS disk. 
iMartin-Zai'di et al] (12008 ') modeled the excitation conditions of 
H2 and the FUSE spectra of HBes stars using the Meudon PDR 
Code and showed that the observed gas was located in relatively 
difiiise regions (10 -3000 cm"-') at significant distances from the 
central stars (0.03 - 1.5 pc). 

None of the four targets with spectra of sufficiently high S/N 
showed any evidence for H2 emission. From the 3cr upper lim- 
its to the emission line flux, we calculated upper limits on the 
column density and mass of H2 for each star. We found that the 
column densities should be lower than ~10^^ cm"^ at 150 K, and 
lower than ~10^' cm"^ at 1000 K, which does not contradict the 
column densities derived from the FUSE observations. Upper 
limits to the masses of warm gas have been estimated to be in 
the range from ~10"^ to ~6 Mjup (1 Mjup ~ lO"-' Mq), assuming 
LTE excitation, and depending on the adopted temperature. 

It should be pointed out that mid-IR H2 lines only probe 
warm gas located in the surface layers of the observed media be- 
cause of the opacity of the interior layers when dust is present. 
The surface layers of any potential disks surrounding our target 
stars do not therefore contain sufficient warm gas to enable it 



to be detected in emission at mid-IR wave lenghs. Following the 
calculations of ICarmona et al. I (120081) . the H2 mid-IR lines pro- 
duced by a gas-rich disk with Tgas = Tdust sh ould not be observ- 
able with existing instruments . As shown bv lBitner et al.l (l2007h 
and Martin-Zaidi et al.] (12007 ^, the warm H2 can be detected in 
the CS disk of a Herbig star, but to explain the detection, par- 
ticular conditions have to be assumed for the gas and dust, such 
as Tgas > Tdusr, which may be created by gas heating by X- 
rays or UV photons. Our non-detections therefore imply either 
that the gas and dust in any potential disks are well mixed and 
have almost equal temperatures, or that these disks are simply 
not present at all. 
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